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Preface 


This final report documents the detailed results of a system application 
study, involving the use of the Hughes Stellar Tracking Attitude Reference 
System (STARS) as the attitude reference for precision pointing (with .001 
degree accuracy) of an EOS spacecraft. The STARS accomplishes precision earth 
pointing by utilizing a single, clustered star tracker assembly mounted on a 
non-orthogonal , two gimbal mechanism, driven so as to unwind satellite orbital 
and orbit precession rates. 

The fundamental objective of the system application study was to explore, 
by means of realistic preliminary design, analysis, and simulation, the limi- 
tations placed upon STARS performance in a typical low altitude, inclined orbit 
earth observation mission by the non-ideal characteristics of the system hardware 
and spacecraft behavior- A basic ground rule of the study was that gyros were 
not to be used to provide short period stabilization. 

The overall conclusion reached in this study is that a relatively straight- 
forward precision pointing system design, utilizing STARS on an EOS type space- 
craft, and considering the characteristics of real sensing and control components, 
c^n achieve the desired 0.001 degree pointing capability. It was found that the 
lower limit on long term pointing precision was determined primarily by the 
accuracy obtainable from the star trackers. Short term pointing precision 
(jitter) was easily reduced to below 0.2 arcseconds by proper control system 
design, and both good jitter stability and excellent transient response character- 
istics were obtained without the use of gyros. Since all analyses and simulations 
so far have shown STARS to be capable of the desired pointing precision and 
stability, the next step in bringing STARS up to the level of a flight qualified 
system is the development of an engineering model and it is therefore recommended 
that this be undertaken next. 
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Introduction 


The requirements of various projected scientific, conmercial, and 
defense oriented satellite missions indicate a strong need for a straight- 
forward, conceptually simple, mechanically and electronically uncomplicated, 
and highly reliable attitude reference and precision pointing system, which 
has the potential of achieving and the possibility of eventually exceeding 
0.001 degree pointing capability. In response to this need, Hughes Aircraft 
Company has evolved the STARS (Stellar Tracking Attitude Reference System) 
concept which avoids many of the problems associated with the precision 
pointing schemes by the use of a compact nniltistar tracker and a simple, 
unconventional two-gini)al assembly. The advantages of this approach are 
as follows: 

1) The stellar reference is a single, clustered star tracker 
asseinbly on a single gimbal mount, thus providing a imi- 
form thermal environment for its components, minimum dis- 
tortion, and involving minimum interaction with other space- 
craft components. 

2) The STARS concept is not dependent upon high precision gyros, 
with their attendant reliability and cost disadvantages. 

3) The star tracker concept provides a high level of redundancy 
in its basic design, being able to track more than the minimum 
required number of stars most of the time. In addition, it pro- 
vides gradual degradation of performance in case of failure of 
one or more individual trackers. 

4) The star tracker design and the particular scanning technique 
used avoid the usual requirements for critical alignment of the 
scanning device to the boresight axis. 

5) The STARS concept operates with essentially constant rates about 
both gimbal axes, thereby avoiding many of the more significant 
sources of angular error normally associated with gimbal systems 
and their drives. 

6) The STARS concept minimizes the amount of computational capability 
required of the on-board computer. 




7) The STARS concept can operate without the conventionally required 
on-board star catalog since the locations of only eight preselected 
stars are required. 

This concept was investigated in NASA funded feasibility study (contract 
- NAS 5-21508) which ended in June 1971. The design configuration, developed 
in the feasibility study, utilized a clustered set of eight star trackers 
mounted on the inner glmbal of a nonorthogonal two-axis gimbal system. As 
illustrated diagrammatically in Figure 1-1, the inner gimbal is held essen- 
tially parallel to the earth's polar axis, while the outer gimbal is parallel 
to the satellite pitch axis. Under ideal conditions, the satellite may rotate 
about its pitch axis while keeping the yaw axis pointed precisely vertical, 
and the star trackers continue to point at their assigned stars while the orbit 
processes about the polar axis. The star trackers are arranged and their signals 
processed so that two are tracking stars at any one time regardless of earth and 
satellite body occultation directions. Physical motion of the star trackers in 
inertial space to permit the tracking function is accomplished by controlling 
vehicle attitude. Precision pointing of the satellite vertical is achieved 
by coimnanding the proper pitch and polar axis gimbal angles, based on ground- 
determined and regularly updated ephemeris data. An on-board digital computer 
keeps track of local satellite time, calculates and sets proper gimbal angles, 
and performs the attitude control signal processing and command functions to 
enable star tracking as required. 

This one year study showed that STARS appears to be a feasible 
attitude reference system and that no theoretical obstacles exist to a 
practical implementation. Systems analysis showed that a set of eight stars 
could indeed be found that would provide a continuous attitude reference for 
any arbitrary orbit. A straightforward attitude acquisition procedure was 
developed, and the on-board computational requirements were established. A 
conceptual design and layout of a suitable gimbal system for the star trackers 
was deVv'^.loped and the availability of the required components was established. 
Servo analysis and simulation of the gimbal positioning systems showed that 
precise positioning at the required position accuracies and rates is theoretically 
feasible. In addition, assembly and test methods for the STARS gimbal system 
were explored as an essential part of an overall feasibility study. 
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A detailed analysis and preliminary design of a star tracker system 
was performed in the feasibility study. This involved a system analysis of 
the proposed tracker system including optics, signal processing, and sunshade 
design. A mechanical design study resulted in a detailed preliminary mechanical 
layout of a tracker cluster and verification that the design concept was feasi- 
ble. In addition a detailed error analysis of the total system, based on antici- 
pated thermal and structural characteristics, showed that the STARS mechanism is 
inherently suitable for high precision attitude reference purposes. 

An area excluded from consideration in the feasibility study was the 
precision pointing capabilityof a real spacecraft employing STARS for sensing, 
since such a study is spacecraft configuration dependent. A preliminary, 
linearized analysis for a typical spacecraft, performed since the completion of 
the feasibility study, has indicated hov^ever that under idealized sensing and 
torquing conditions, with pure angular rate feedback available, precision verticd 
pointing to 0.001 degree is practical. The questions that were not answered 
at that time, and that required further detailed study, concerned the limitations 
placed on precision pointing by realistic star sensor characteristics, signal 
processing requirements, control system bandwidth limitations due to sensor 
noise, and the requirement to generate attitude rate damping information from 
star sensor signals rather than a built-in gyro. These are the principal 
questions that were addressed by this system application study. 

The study can be divided essentially into three parts: star tracker investi- 

gations, rttuftude control system design, and performance analysis and simulation. 
Each of these areas is discussed in turn on the following pages. 
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2. Star Tracker Performance 

In order to develop more confidence in the ability of the STARS tracker 
to achieve the required performance in the postulated environment, the tracker 
noise model developed during the feasibility study was re>examined. Concurrently, 
to provide improved insight into photomultiplier noise characteristics obtained 
experimentally on another program, and to confirm the analytical approach, ad- 
ditional tests were performed on three phot omul ftp Her tubes. The analysis and 
che experimental results are discussed below. 

2.1 Analysis 

The star detection process basically consists of setting a threshold level 
and defining voltages below that level "noise" and voltages above that level as 
"signal". The determination of the appropriate threshold level setting requires 
a statistical evaluation of the probability of a noise pulse being greater than 
this lev'el (false alarm) and a signal + noise pulsrj being smaller them this level 
(missed pulse). This analysis ic based cn the assumption that a false alarm is 
equivalent to a missed pulse in respect to tho perturbation it causes to the 
system's tracking loop. 

The STARS star tracker uses photomultiplier tubes as optical detectors. 

The model that best describes the pulse amplitude distribution of the photomulti- 
plier tube's noise output has been the principal subject of this analysis. It is 
this model that permits a determination of where to set the detection threshold and 
prediction of the false alarm and missed pulse rates for this threshold. In earlier 
analyses, gaussian statistics were used to describe both the signal and noise 
characteristics of the PMT output. Experience on the Space Precision Attitude 
Control System (SPACS) star sensor program found that this assumption yielded 
extremely optimistic results in comparison to laboratory measured performance. 

The noise output of a PMT is primarily a result of PMT leakage current, star 
background irradiance, and illumination from off-axis bright sources such as sun, 
earth or moon being reflected and scattered by the sun shade into the field of 
view. Of these, the reflected and scattered solar illumination was considered 
to be the dominant performance limiting noise source as discussed in the STARS 
final report. 

I addition to these sources of noise; space radiation, in particular high 
energy galactic particles, may contribute to the PMT > iitput noise. 
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A study that was perforcied to explain the discrepancy between the 
measured and expected false alarm rates on the SPACS program found that the 
amplitude of the PMT output noise pulses exhibited an exponential distribution 
rather than gausslan. This distribution results in much higher false alarm 
;uates as will be shown. 


The present analysis Includes an examination of the pulse amplitude 
distribution of three representative photomultiplier tubes operated at the 
conditions of background Illumination and electrical bandwidth anticipated 
for the STARS sensor. 


The solar irradlance that is scattered and reflected from the sun shade 
will result in an average (d.c.) output at the PMT anode. This current can 

D 

be related to an equivalent signal irradlance at the entrance aperture of the 
optics by 




1 

R G A 


B 

eff 


where 


R = Photocathode Responsivity 


G = Photomultiplier Gain 
A^^^ = Effective Optical Aperture 

The variation of background noise pulse amplitudes is generally described 
by a gaussian distribution for which the Iq standard deviation or rms noise is 
given by 

ijj(la) = J2e ^ ^eff ^ 


where 

6 ® average gain of an individual dynode 
G “ multiplier overall gain 
4 f - noise equivalent bandwidth 


The effective signal irradlance equivalent of this noise current is giver. 


by 


n 


A ,, RG 
eff 


-2 

(watt cm ) 
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The baseline design parameters, applicable to this analysis, for the 
STARS sensor are shown In Table 2-1. 


TABLE 2-1 BASELINE DESIGN 
Photocathode Responslvity, R =» .075 amp/watt 

2 

Effective Collecting Aperture, “ 6.4 cm 

Noise Equivalent Bandwidth, ^f = 350 hz 
Gain of a single dynode, 6 = 2.5 

-13 2 

Background Irradlance, * 1.1 x 10 w/cm 

-13 2 

Minimum STAR Irradlance, = 1*1 x 10 w/cm 

Pulse Visibility Factor, V^/Vas = .75 


Normalizing Irradlance levels and noise currents to the appropriate value 
of CT provides a simplified means of comparing data taken under different con- 
ditions and relating It to the analysis using baseline design parameters. 

For the above design parameters, it is found that ^ = 6.52 x 10 watt 

-2 

cm . The minimum star therefore can be expressed as 


^eff ■ “ 12.65(7 


For a gaussian distribution of noise pulse amplitudes, the probability 
of a pulse exceeding the mean value by an amount T can be described by 


P^^(T) 


1 _ 



exp 



dT 


The gaussian function is generally used to describe noise amplitudes greater 
than and less than the mean value with equal probability. For the purpose of the 
analysis of false alarms, only positive amplitudes need by considered, resulting 
in the j factor at the beginning of the above equation. 

If the output noise exhibits exponential behavior the distribution of noise 
amplitudes greater than the mean level can be described by 
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A(T) * ^ exp ^ 

If the assumption is made that the output noise amplitude is greater than 
the mean level half the time, then the probability of a false alarm as a 
function of threshold setting above the mean level can be expressed by 



In the baseline STARS concept, each photomultiplier tube is only active 
during two 90° windows. This results in a factor of 2 reduction in the 
number of effective false alarms. 

The probability of missing a pulse is described by the probability that 
the amplitude of the signal plus noise will be less than the threshold level. 
One must be certain to include the noise because of the finite probability 
that the direction of the noise pulse at the time of a signal pulse will sub- 
tract from the signal thus increasing its probability of falling below thres- 
hold. Since both positive and negative deviations of signal and noise from 
their mean values must be considered in this portion of the analysis, use of 
a gaussian distribution is most appropriate. Data obtained on the SPACS 
system supports use of the gaussian model. 

The standard deviation or rms value of the signal + noise is given by 

‘n ' Kit ' <“.tf "b> 

and the corresponding signal irradiance by 

^n -2 

c' “ ^ (watt cm ) 

eff 

For the baseline case in which = ,/Tg 

The missed pulse rate is given by 
MPR = Pjjp . N 

where N is the normally expected rate of signal pulses which for a STARS 
sensor is four per second. 
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2.2 Experimental Results 

A test was set up In which a light source provided a background irradiance 
at the PHI photocath.ide. The d.c. or average and the rms noise currents were 
measured at the output of a 300 hz (3 db) cutoff filter. A threshold was set 
at the average and the number of frequency of pulses exceeding this value were 
measured on a counter. This zero threshold (T * 0) level was found to yield 
a false alarm rate that could be approximately described by 

FAR as Af 

The earlier SPACS study reported a similar finding^ however, certain as- 
pects of that study appear to indicate that noise amplitudes below the mean 
noise level were neglected, which may imply false alarm rates a factor of two 
greater than those observed in these tests. 

The threshold voltage level was increased and data was collected of the 
frequency of pulses that exceeded each threshold level. These frequencies or 
count rates were divided by the zero threshold rate to provide the probability 
of false alarm as a function of threshold and normalized to a value of .5 for 
the P^^ (T = 0). 

The three tubes used in this test had the sune type of photocathode and 
were manufactured by the same manufacturer. One tube, 7762E, has a slightly 
larger photocathode area than the other tubes and was manufactured around nine 
years ago. A second tube, 13873, was manufactured around four years ago and 
has been used quite extensively for SPACS engineering tests. The third tube, 
19391, is of very recent manufacture and has not been used other than for 
acceptance type tests. 

As seen in Figure 2-1, one tube (13873) behaves in an almost perfect expo- 
nential manner and the other two tubes are closer to gaussian. P^f^ 19391 shows 
a tendency toward exponential behavior for large amplitude pulses. 

The nature of the measurement techniques used to obtain the above data 
necessitated using data taken over extended periods, many days in some cases. 

This caused some concern over the accuracy of applying the data to short time 
intervals. A solution to this shortcoming was found to be through the use of 
a multichannel analyzer that samples noise and counts the distribution of samples 
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as a function of their an^litude. 

A Northern Scientific, Model NS'900>1A multichannel analyzer was employed 
for this direct photomulitplier dark and background distribution analysis. 

In addition a Tektronix 347 oscilloscope with a 1A7A plugin was utilized to 
pruuuce a 4B coincidence gate allowing the analyzer to be operated in a multi- 
channel A/D conversion mode for analysis of the slowly varying filtered wave- 
forms. 

An examination of the data collected in this test (See Figure 2-2) con- 
firms the previous conclusion that the distribution of large amplitude pulses 
dif fes markedly between tubes of the same family. Even after substantial fil- 
tering (300 hz bandwidth) some of the photomultipliers show an exponential 
form of tail distribution that is noticeably affecting the background distribution. 

The exact cause of this variation in behavior £imong tubes is not known at 
this tine. It is speculated that the presence of gas in the tubes may significantly 
affect these distributions, in particular, the large amplitude pulses. 
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2.3 Results and Conclusions 


Figure 2>3 describes the several probability function distributions 
that have been discussed In this report for both the noise and signal + 
noise cases. Using these curves, the false alarm and missed pulse rates 
were confuted for the exponential distribution as a worst case condition 
and for PMT 19391 as a more optimistic case that can probably be realized 
through photomultiplier tube selection. Tables 2-2 and 2-3 sunanarize 
the predicted false alarm and missed pulse rates as a function of threshold 
setting for the above two cases. 

The signal to noise ratio of the STARS sensor as is used In the deter- 
mination of its sensitivity and stability is given by 


SNR 


“.ff VV 


For the baseline design parameters, a SNR of 8.95 is obtained for the minimum 
required star. 

It must be kept in mind that the performance characteristics described 
are a function of the angh between the pointing direction of the star tracker 
and the location of bright sources, which In this analysis was 15° for the 
earth and 30° for the sun. Performance %n.ll improve at greater angles until 
other sources become the limiting noise factor. In the SPACS program it was 
determined that galactic particles probably provide this lower limit. 

In the course of the experimental portion of this study, the background 
source was turned off and a count was made of false alarms from all other 
sources. For a typical threshold setting, the count rate was around two orders 
of magnitude lower than that obtained with a bright source background. 

In summary, this study indicates that as a worst case situation, a 
false alarm and missed pulsed rate of 75 per hour can be anticipated. A more 
optimistic, yet realistic condition predicts around 4 false alarms and missed 
pulses per hour. It should be noted from the probability values listed in 
Tables 2-2 and 2-3 that the occurrence of two consecutive perturbations has 
very low probability. It is also very obvious that an order of magnitude 
improvement is obtainable if desired (at the expense of electronic complexity) 
if time gating of the star pulses is applied as part of the signal processing. 
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FIGURE 2-3. PROBABILITY OF FALSE ALARM AND MISSED PULSE AS FUNCTION 
OF THRESHOLD LEVEL 
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TABLE 2-2 


False Alarm & Missed Pulsed Predictions. Exponential 

Case 


Threshold 

^’fa 



^MP 

MPR, hr"^ 

7a 

4.5 X 10 ^ 

565 


3.5 X 10"^ 

.505 





-4 


8a 

1.65 X 10 

208 


5.8 X 10 

8.35 


-5 



-3 


9a 

6 X 10 

75 


5 X 10 

72 


-5 



-3 


9.2a 

5 X 10 

63 


7.2 X 10 

104 


-5 



-2 


9.5a 

3.5 X 10 

44 


1.3 X 10 

187 


TABLE 2-3 

False Alarm & Missed Pulse Rates, PUT 19391 


Threshold 

"fa 

FAR, hr"*^ 

— 

"mp 

MPR, hr-1 

6a 

2.4 X lO"^ 

30 


1.0 X 10"^ 

.026 


-6 



-5 


7a 

6.6 X 10 

8.2 


3.5 X 10 

.505 


-6 



-4 


7.5 a 

3.5 X 10 

4.4 


1.5 X 10 ^ 

2.16 


-6 



-4 


7.7 a 

3 X 10 ° 

3.8 


2.2 X 10 

3.2 


-6 



-4 


7.8a 

2.75 X 10 

3.5 


3.0 X 10 

4.3 


-6 



-4 


8a 

2.2 X 10 

2.8 


5.8 X 10 

8.35 
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3. Attitude Control System Design 
3.1 Spacecraft Dynamics 

The spacecraft dynamics used In the simulation Included nonlinear 
dynamics to obtain a more accurate representation of the spacecraft dynamics 
for the STARS application. The hl^ degree of accuracy desired In the 
dynamics model complements the high degree of accuracy expected from STARS. 
Nonlinear terms as well as product of Inertias are Included In the analysis 
to formulate a more representative dynamics model of the spacecraft. 

The Inertia matrix for the spacecraft Is shown below. The coordinates, 
X, y, z are inertially referenced general orthogonal coordinates. 



where 

^ij =-J i 

^ii J (j^ + 

1» j» = coordinates, origin c.g. 
m = mcds 

The angular velocity and acceleration matrices for the vehicle are 
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FIGURE 3-1. BLOCK DIAGRAM FOR VEHICLE DYNAMICS 
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where 6 = angle rotation about q-axls, q = x, y, z 

The moment torque equation is 

• 

= H + u X H 
eg ” 

where H = I ^ 

The moment equation for the EOS vehicle can be calculated to be (See Appendix 5-1) 
• • 

^ ‘“r ^RY ^RY '“P “R + ^Y '“y ‘ <»y 

^ “P ■ ^RY <"1 ■ ^Y “'r ^R “r «*Y ^RY '4 

^ ^RY ^ ^P “'R “>P ■ ^R “»R <"P ■ ^RY '"P ^ 

Figure 3-1 shows the block diagram for the spacecraft dynamics based on 
the set of equations above. The torque exerted on the spacecraft produces 
roll, yaw, and pitch angles 0, Y, and 0, respectively. The torque input is 
supplied by the control system that receives attitude error information 
from STARS. The spacecraft attitude angles, 0 Y, and 9 are then reflected 
in the errors sensed by the star sensor/telescope. 

3.2 Disturbance Torques 

Sources of disturbance torques on an earth s' tellite are environmental 
(external) or internal. Those sources, due to interactions of the spacecraft 
with the environment, are gravitational gradient, solar-radiation pressure, 
aerodynamic pressure, meteorite impact and magnetic interaction between the 
earth's field and spacecraft electric currents. Internal sources are the 
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expulsion of aatter (e.g., thrusters), activation of spacecraft elements 
or appendages and Interaxis coupling. 

Disturbance torques can be either cyclic or secular as referenced from 
inertial space. Moc ntum storage devices (e.g., reaction wheels) can be 
utilized to accommodate angular momentum due to cyclic disturbance without 
requiring unloading and exceeding minimal attitude error. Secular torques 
require unloadi g of stored momentum periodically because of continuous' 
momentum buildup. 

SOLAR PRESSURE TORQUE 

The resultant solar pressure torque is a function of the effective 
area of the spacecraft exposed to the sun, the reflectivities of the 
exterior surfaces, and the distance between the center of mass and center 
of pressure. 

The basic equations for solar pressure forces are; 



- 


f 

A A 

/ (1 + VjKj) 

V • n 


(y 


n) n d A. 


A. 

1 


(Normal Forces) 


+ 


f 

A A 

1 

V * n 


(n X V X n) d A 


I 


(Shear Forces) 


+ 



f “1“ 



A 

V * n 


n d Aj 


(Diffuse Reflection 
Forces) 
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where the terms are given as; 

= reflectance coefficient over the ith surface 

= diffusivity coefficient over the surface = 1 for specular 
reflection) 

n = unit vector directed outward normal to dA^ 

= solar radiation pressure constant for the normal incidence (v.ft) 

- 1 and complete absorptivity (v = 0) = 9.A x 10~® pound/foot^ 

= the area of the ith surface 

^ = unit vector directed to sunline 

The torque produced by the solar radiation pressure is 

Ts ^ 

where is a vector from the vehicle center of mass to the surface (dA^) 
center of pressure. 

Normal force effects were considered adequate for the generalised EOS 
spacecraft configuration used in the study. The solar panel is sun oriented 
and contributes heavily to the normal force, while the remaining surfaces tend 
to produce cancellation of shear force and/or have low effective reflectances. 

The maximum solar pressure torque has been calculated to be 1.0 x 10 ^ 
foot-pound for the 9:00 A.M. orbit. 

For the computer simulation, the sun line is assumed to be along the 
vernal equinox and the spacecraft is assumed to be symmetric about the pitch 
axis such that the center of pressure is located along the pitch axis. Torque 
about each axis is approximated by a sinusoid based on the orbit rate. The 
approximations used are: 


?A 



\oll " 

Tpi^ch * foot-pound 

T„ “ 1.0 X 10 ^ sin <u I 
Yaw 0 


t foot-pound 


foot-pound 
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where 


_3 

ID^ » 10 radlan/second 

The solar torque presented above should provide a representative 
disturbance torque for observing STARS perforaance. 

AER(»YNAMIC PRESSURE TORQUE 

The aerodynaaic pressure effect on the spacecraft at an altitude of 
1000 ka produces a torque on the order of 2 x 10 ^ foot-pound, priaarlly 
in the pitch direction for an angle of attack of eight (8) degrees. The 
general equation for aerodynaoic torque is 

^ *A <*cg ■ %) » 


where 


= torque due to aerodynamic pressure, foot-pound 


= atsiospheric pressure, pound/ foot 

2 

= surface area exposed to aerodynaadc pressure, foot 

a = vdiicle angle of attack, radians 
Compared to solar pressure torque, aerodynamic pressure torque is small. 


GRAVITY GRADIENT TORQUE 


A spacecraft with non-equal principal momenta of inertia tends to 
align itself in the gravity field so that its principal axis of minimum 
inertia lies along the gravitationa] vertical. The torque components due 
to gravity gradient are 


^g roll 


Tg pitch 


T yaw 
8 


^ »o “p - *Y> ^ * 

i -o - S’ ^ * 
0 
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%fhere 9 ? 10 ^ radian/second orbital rate 

Ir, Ip, are the v^icle oioaients of inertia about roll, pitch and yaw, 

2 

£oot-pound>aecond ; 0 and 0 are the roll and pitch error angles, radians. 

The gravity gradient effect produces a torque about the pitch axis 
-8 

equal to 2.6 x 10 foot-pound, which is small compared to the solar 
pressure torque. 


MAGNETIC PRESSURE T(«QUE 

The magnetic field of the earth exerts magnetic pressure on the space- 
craft through an interaction with electric currents on the spacecraft. The 
torque applied to the spacecraft is given by 

T = 6.86 X 10 ^ M X B foot-pound 

where 

— 2 
M = magnetic moment, amp- turn- foot 

B = magnetic flux density due to earth's magnetic field, gauss 
The magnitude of the earth's magnetic field can be expressed as 


where 


B = 


1.43 X 10 


10 


(sin Y -t- 4 cos y) 


1/2 


r = distance from dipole center to field point, radius of orbit 
(nautical miles) 

Y = angle r makes with the earth's north-south dipole centerline 

(radians) or « = m t where w is the orbit rate. 

* o 

The distance r is the sum of the earth's radius (3444 nm) and h^ the 

altitude of the orbit. For an altitude of 540 nm, the magnetic flux 

-3 2 

density of the earth is 4.6 x 10 gauss. Assuming 600 ampere-turns-foot 

for the spacecraft's magnetic moment, the peak value of magnetic pressure 

torque is 1.9 x 10 ^ foot-pound. The torque varies with position along the 


orbit. 
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3.3 Coapematlon 

The vdilcle control subsTStea that receives attitude reference data 
froai the STARS senaor and giabal subsysteas aust provide the proper 
coapensation for the overall STARS attitude control syster' in order for 
the vehicle's response to be stable and have the desired pointing 
characteristics. The coapensation that provides the control law for this 
system has two major tasks; it aust hold the vehicle pointing error due to 
disturbance torques to within +.20 arcseconds, and it aust also attenuate 
the STARS sensor noise. The disturbance torques acting on the vehicle 
occur at orbit frequency, tdiile the STARS sensor noise occurs at 
frequencies a couple of magnitudes hi^er than orbit frequency. 
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A simplified block diagram of the overall STARS attitude control 


system is shown below. 



The requirement that the vehicle pointing error due to maximum 
disturbance torques is less than .20 arcseconds; places a magnitude 
contraint on the transfer function, 9/T 


Q (1 tti) 


Tn (j s.) 


10 


-6 


10 


N 

1 radian 
10 foot-pound 


From the figure 

9 (j a») 




i„ (j «) 
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1 + Gj^ (j (u) (j oi) 
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This expression yields a gain of 10 for at orbit frequency. 

To reduce pointing error due to STARS sensor noise » the gain of 0/T^ 
should be made as small as possible at the higher frequencies. With 
these gain constraints and other considerations in mind, the following 
compensation was chosen: 


( 8 ) 


K (T^ 8 + 1)^ 
8 (t^ S + 1) 


K = .01 foot-pounds /radian 
= 67 seconds 
= 6.7 seconds 


The gain K is chosen, as described above, to achieve the required 
pointing accuracy in the presence of disturbance torques. The double 
lead is required to obtain sufficient gain and phase margins for system 
stability. An integrator was added to alleviate a steady-state error, 
a result of the annual precession of the orbit plane about the center of 
the earth, which was encountered during simulation. The lag is used to 
produce a roll off at higher frequencies to filter out STARS sensor 
noise. 


3.4 Control System Configuration 

The control system torques the spacecraft to produce attitude 
changes. Torquing can be produced by momentum exchange between the 
vehicle and a momentum storage device, such as the reaction wheel. Three 
reaction wheels can be used to handle three attitude angles of the 
spacecraft. 
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TORQUE 



FIGURE 3-2. CONTROL SYSTEM CONFIGURATION 
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Other means of applying torque are available, but for the sample 
example of torqulng separately In each of the three axis and of making 
use of the cyclic nature of torque disturbance from the environment, 
the reaction wheel torqulng system is a good candidate. Jet thrusters 
are also used, and they transfer momentum from the wheels during the 
unloading of the wheel momentum. 

The control that produces signals to command the reaction wheel 
Includes a shaping network or compensation to give good control 
characteristics to the control systm. The compensation is designed to 
attenuate sensor noise and to enable the control system to maintain 
+.20 arcsecond vehicle pointing accuracy under the effects of disturbance 
torques . 

The simulation uses this relatively simple control law. The 
results show that this control law is effective in providing good 
response characteristics in ontrolllng the nonlinear spacecraft dynamics. 
Figure 3-2 shows the control system configur '.tion. Included in the figure 
is the logic box for unloading the wheels. Various torques that the 
vehicle experiences are also shown. 

The use of a relatively simple control law was desirable to show 
that STARS can operate with a simple control system rather than only 
with a limited and/or complex type of control law. The use of a simple 
control law also allows for the use of more specific or precise control. 

If needed. A more complex control law can be designed to handle more 
difficult control problems If such problems should arise. 
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4 


Performance Analysis and Simulation 


4.1 Preclalon Gtmbal System 

Figure 4-1 shows the configuration of the precision glmbal pointing 
control system. This configuration is Identical for both the pitch and 
polar axes. 

The command to the system (9^) is a binary number that reprerents 
the desired gimbal position which is calculated on board and updated at 
predetermined Intervals using ephermis data. The actual position of the 
gimbal is monitored by use of an Inductosyn. The output of the Inductosyn, 
a binary signal, is fedback and a binary subtraction is made between this 
signal and the conmand. The resulting error signal is passed first through 
a digital to analog converter and then through a shaping network. The 
output of the shaping network is a voltage that drives an Inland brush dc 
motor which positions the gimbal shaft. The positioning of the gimbal 
shaft by the motor is corrupted by bearing noise. This bearing noise which 
is due to mechanical tolerances and some uncertainties is described in 
Table 4-1. The other noise source in the gimbal system is on the Inductosyn 
measurement of the actual gimbal position. This is described in Table 4-2. 
Both noise sources are assumed to be white with zero mean. 

4.1.1 Pitch Gimbal Simulation 

The pitch gimbal servo is used to unwind the orbital rate of the space 
craft about the earth; this must be performed in order to maintain the pitch 
axis of the telescope fixed in Inertial space. The conmiand angle to the 
pitch servo is the orbital rate of the spacecraft about the earth. This 
command is approximately 240 arcseconds per second and is updated every .015 
seconds. Due to this relatively high angular rate and to the high load 
inertia of the pitch axis, a stick-slip friction model can be avoided and 
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FIGURE 4-1. PRECISION GIMBAL SYSTEM 
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TABLE 4-1 
BEARING NOISE 


— 

Ittm 

1 

F.rfor 

CUBt* 

Error Qiunllty 

Definition 

Aaalgned 
Value 
Maximum, 
mic roinches 

Angular** 
Contribution 
in Aaaembly, 
micro radians 

1 

Rema rfca 

1 

E 

RE 

L URB-lROe 
t (IRB.CRC)* 

Eccentricity. Inner race 
'groove u> inner race bore 
fTIRl 

Rceiduel error of 1 

bO 

IS 

Not 

applicable 

Each inner race marked with high spot. 
High spots of bearing psir aligned, thus 
causing residual crrui due tu a) ntagni^ude 
difference between race eccentricities and 
b) deviation in true high spot location and 
alignment. 


U 

L (IRB.lRCle 

Uncertainty of 1 

15 

T5‘ ' * 

Uncertainty in (E) i- URB-lRG)c 

1 

U 

L (IRRF-IRC)AW 

Orthogonal component (radial 
direction). Axial wobble of 
inner race groove to inner 
race reference face 

7.5 

10 

Axial wobble high spot rarely coincides 
with eccentricity high spot so that cor- 
rection cannot be n.adc simultaneously 
for wobble. Assumed worst case for 
opposite bearing pairs dismetrically 
opposed; howrvi.r. mvan was asiunied 
within each boring pair. 

S 

V 

L (B-B)V 

• 

Variance. BalUto-ball 
diametcra (for bearing 
ball complement) 

5 

O 

o 

212 

Assuming: larger balls (!/2 ui the 
complement) dianTeiricxlly oppose sn^aller 
balls. Worst possible case bath boa ’’ing 

pairs. • 

■ 

u 

L PC 

Sisc of particulate 
contamination 

95 

15. 1 5 
10 

Assuming: some particles at tirrics enter 
ball to race contact and not necessarily 
aimultancously in both bearings and pairs 
in diametrically opposed locations. 

■ 

E 

L (Stf.IRK 

Fit between inner race 
bore and shaft 

Inter- 

ference 

”1 

Interference to exist over entire opera- 
tional thermal range. 

H 

E 

L (H-OR|f 

Fit between outer race 
and housing 

Inier- 

ferencr 



1 

E 


Shaft eccentricity between 
bearing locations 1 and 2 
(TIR) 

60 

Not 

applicable 

Combined with item 1, rucet are S'>.;'aUrlv 
positioned to minimirr^ effective to’^L 
eccentricity; i.e.', shaft raiational axis 
defined by race grooves. 

10 

RE 

L SH 

Residual error of shaft axi* 
definition (1 and 9 
combination) 

10 

10 * 


11 

V 

L r.SH 

Uncertainty of shaft axis 
definition (1 and 9 
combination) 

15 

15-- I 5 
10 ’ 


12 

U 

t «M-v9)) e.SH. 

TH 

Uncertainty in 1 and 9 com- 
bination due lo ronho’i\o- 
geneous response to thermal 
cha nges’^<='= 

2.5 


Thermal diffcrentisl exoangion 
change in initial i nto rfe rvnee ftif: u ay 
not take place equally evrrvuhere due tu 
nonhomogencity of structurts and matcrtals. 

U 

u 

TH 

Balt contact position change 
due to outer race to inner 
race Ihcrn^al condition 
change 

2.5 

^•D.SO 

Even if both inner and outer races experience 
identical change in temperature (no gradient) 
and if both bearings experienced identical 
change* larger outer race will change at a 
higher cate than inner race. Contact an>{le 
shifts. Assume; bearing pair (1) experiences 
thermal shift larger or smaller than that 
experienced by bearing pair f!) this causes 
angle uncertainty (11) due to difference in race 
curvatures and sites between btarjn^ 

pairs (1) and (2). 

U 

u 

L (OR-H:$G}«, 
TH 

Housing to outer race 
Induced eccentricity change 
due to nonhomogencous 
interference fit change due 
to thermal variations 

2.5 

A. 0.50 

Original outer race ana housing errura (as 
aaaembled) arc tumped in an offset an<le 
error — a constant. However, if changes 
occur, the otfse* an^lr changes in an 
unpredictable manner. 

IS 

1 

u 



Ips'LE) 

Uncertainty due to launch 
environment ceused Brincll 
• pota 

30 


75% Brinell apoi depth of indent.-ition in one 
bearing, taking total load in radial direction, 
(60 g peak acting for > 1 second). 

error* RR: rcfidu«l error, end U; uncerUinty. 

^-1 microradian * 0. 20^arc»ec. 

;£jo!k temp;: rit'i re N’arialion Unntc«' Ij ♦ lO^F; radial vraditnr 

to 5^i‘ 
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TABLE 4-1 
(CONTINUED) 


A 

Item 

Error Quantity, 
micro radians 

■ 

Quantity 
Squared, 
micro radians 

Extension 



summary 


2 

1.5 

1.5 

2=5 (errors) 

10 

3.0 

9.0 

Rss = 3.24 microradians 




~ 0. 7 a rcsec 

3 

1.5 

2.25 

2= 37,8 (uncertainties) 

4 

• 1. 5 

2.25 



Rms = /37. 8 

5 

1.00 

1.00 

V Q-1 

6 

4. 5 

20.3 

2 = 2. 17 microradians 




~ 0. 45 a rcsec 

11 

1. 5 

2.25 


12 

0. 50 

0.25 


13 

0. 50 

. 0.25 


14 

0. 50 

0.25 


15 

3.0 

9.00 


Therefore: Total error and uncertainty contributions (maximum) = 

- 

0. 70 + 0. 45 

= 1.15 arcsec 

• 
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TABLE 4-2 


INDXTGSYN ERRORS 


A. 

Errors systematic with poles 
(512 per revolution) 

0.45 arcseconds 

B. 

Errors systems per vevolution 



1. 

Centering indication on Inductosyn discs 

0.68 


2. 

Effect of bearing bias 

0. 50 


3. 

Rotor to housing alignment 

1.00 


4. 

Stator to shaft alignment 

1.00 


5. 

Wobble 

0. 20 


6. 

Bearing random effect 

0. 10 

RSS TO TAL 

1.72 arcseconds 
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a conventional coulomb friction model used. 


A separate computer simulation of the pitch gimbal servo was per- 
formed using the MIMIC language. The computer listing of this pitch 
simulation can be found in Appendix 5.2. Because the command angle to 
the pitch gimbal is updated every .015 seconds, a very small integration, 
step is required, thus only short simulation runs (approximately 20 
seconds) were economically feasible. The pitch ginobal system response is 
shown in Figures 4-2, 4-3, and 4-4. In all three figures, the independent 
axis is scaled to seconds of time and the dependent axis is scaled to arc- 
seconds . 

Figure 4-2 shows both the command angle (.dashed line) and the actual 
angle (solid line) as a function of time. Both curves lie on top of one 
another, making the gimbal error barely visible on this scale. Figure 4-3 
shows the difference between the command angle and the actual angle on a 
much smaller scale. The gimbal error is seen to stay almost completely 
within a plus or minus two arcsecond error band. The RMS value of the 
gimbal pointing error is shown by Figure 4-4 to be 1.2 arcseconds. This 
pointing error is obviously well within the 1.80 arcescond error budget 
allotted to each gimbal system. 

4.1.2 Polar Gimbal Simulation 

The polar axis servo is used to unwind the annual precession of the 
orbit plane, which must be performed in order to maintain the polar axis 
of the telescope fixed in inertial space. The command to the polar axis 
servo corresponds to a one revolution or 360 degrees per year rate. This 
command is updated only every 90 seconds. Unlike the pitch gimbal, the 
slow rotation and small angle incremental motion of the polar gimbal cause 
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an elastic behavior between the moving surfaces of the glmbal. As a 

consequence, the polar glmbal slotulatlon utilizes an elastic friction 

model, based on experimental data obtained on another Hughes program. Figure 4>5 

shows the elastic friction characteristics observed on the Hughes OSO 

Program. The STARS polar glmbal operating region Is shown In the lower 

left hand portion of the figure. For this region, the friction torque is 

directly proportional to the angular displacement. The elastic friction 

model that was used In the sla.ulatlon Is shown In the right half portion of 

the same figure. The elastic friction torque in the simulation model is 

equal to K, times A9, when a9 Is not equal to zero, and zero when A3 is 

equal to zero. Where A3 is the Incremental change in the command at each 

update interval. 

A separate computer simulation of the polar glmbal ser\’o was per- 
formed using the MIMIC language. The computer listing can be found in 
Appendix 5.2. Figures 4-6, 4-7, and 4-8 show the polar gimbal system 
response. In all three figures, the independent axis is scaled to seconds of time 
and the dependent axis is scaled to arcseconds. 

Figure 4-6 shows both the command angle (dashed line) and the actual 
angle (solid line) as a function of time. The gimbal angle is seen to 
follow the staircase type motion of the command angle; however, its ampli- 
tude is corrupted by noise. Figure 4-7 shows the error between the 
commanded angle and the gimbal' s actual position. Most of the error is 
seen to fall into a plus or minus 2 arcsecond region about zero. Figure 
4-8 shows the RMS value of the pointing error to be 1.80 arcseconds, which 
is exactly the allotted amount of error for the polar gimbal system. 
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The pitch and polar ginbal slaulations deacrlbed above deannstrate 
that the preciaion ginbal of Figure 4-1 can aeet the hlgli accuracy require- 
ments of STARS. 
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FIGURE 4-2. PITCH GIMBAL RESPONSE 
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FIGURE 4-3. PITCH GIMBAL ERROR 
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RMSriTa. GIMIAL ERROR. ARCSECONOS 



FIGURE 4-4. RMS PITCH GIMBAL ERROR 
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ANGULAR DISPLACEMENT, ARCSECONOS 


FIGURE 4-5. ELASTIC FRICTION CHARACTERISTICS 
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POLAR GIMBAL RESPONSE, ARCSECONHS 
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FIGURE 4-6. POLAR GIMBAL RESPONSE 
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FIGURE 4-7. POLAR GIMBAL ERROR 
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RMS POLAR GIM8AL ERROR, ARCSECONDS 



FIGURE 4-8. RMS POLAR GIMBAL ERROR 
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4.2 Control System 


Figure 4-9 shows the control system configuration b '.ock diagram. 

The command to the control system is the desired spacecraft attitude. The space- 
craft pointing error, which is obtained by comparing the reference attitude 
to the actual spacecraft attitude measured by STARS, is passed through a 
shaping network that commands a torquer to torque the spacecraft in 
order to produce the required attitude correction. This system contains 
two major disturbance sources; the spacecraft disturbance torques and the 
sensor noise on the STARS measurement of spacecraft attitude. The shaping 
network is designed to attenuate the STARS sensor noise and compensate for 
the disturbance torques. 

The control system of Figure 4-9 was simulated under the assumptions 
that the torquer is ideal (no dynamics) and that STARS gives a perfect 
measurement of the spacecraft attitude. The purpose of this simulation was 
to investigate the system's ability to control the spacecraft dynamics and 
compensate for the disturbance torques. 

Figures 4-10 thru 4-12 show the control system response when no distur- 
bances are acting on the system. Figure 4-10 shows the spacecraft 
attitude o”er one orbit to be perfect. Figures 4-11 and 4-12 show the 
spacecraft's body rates. From the definition of the coordinate systems 


(see Figure 4-13), the body rates for perfect pointing should be: 


01 =0slnY8ino^ = 

yaw 


.04 sin O' t 


arcsecond 

second 


= 0 sin V cos o = .04 cos a t 


arcsecond 

second 
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From examination of the three figures, it can be stated that under 
ideal conditions and no disturbances, the control system can point the 
spacecraft perfectly. 

Figures 4-14 thru 4-16 show the control system response when distur- 
bance torques are acting on the system. Figure 4-14 shows the spacecraft 
attitude time history for one orbit period. The large peak in the roll 
response is an initial transient due to the manner in which the simulation 
was Initialized. The pointing error caused by the disturbance torque is 
well within the prescribed accuracy of STARS. Figures 4-15 and 4-16 show 
the spacecraft body rates. The body rates are seen to be unaffected by 
the disturbance torques. 

It can be concluded from the above simulations that the control system 
design of Figure 4-9 can control the spacecraft dynamics and also compensate 
for any disturbance torques acting on the spacecraft . 
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FIGURE 4-9. CONTROL SYSTEM CONFIGURATION 
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0. 0.65 1.8 1.95 2.6 3.25 3.9 AM 52 S.85 6.5 

TIME, THOUSANDS OF SECONDS 

FIGURE 4-10. SPACECRAFT ATTITUDE, ZERO DISTURBANCE TORQUES 
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FIGURE 4-13. STARS GEOMETRY 
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0 0.65 1.30 1.65 2.6 3J5 i.9 4.55 5.2 SJ5 9.9 

TIME. THOUSANDS Of SECONDS 

FIGURE 4-14. SPACECRAFT ATTITUDE. DISTURBANCE TORCKJES PRESENT 
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4.3 Star Tracker 


A block diagram of the tracker model illustrating the relationships of 
its component functions is shown in Figure 4-17. 

The time interval (At) bet%ieen a position reference pulse coming from a 

pickoff on the optical wedge and the pulse produced by the star image crossing 

the recticle slit within the telescope is stored in a register for use by the 
tracker system. The register is updated every 250 milliseconds. A sinusoidal 
relation exists between the time interval and the corresponding error angle, 
given the scan circle radius and its period. Thus, the star line-of-sight 

error angle is derived directly from the pulse time interval measurement. 

The realistic tracker model includes inputs for the two major disturbances 
to the system; missed (or false) pulses and pulse jitter. Both photon noise from 
the star bacl;ground and dark current noise of the photomultiplier tube will cause 
missed (or false) pulses. This subject was treated in Section 2 of this report. 
With respect to this type of disturbance, the following assumptions were made in 
conjunction with the simulation model: 

1) Missed or false pulses will perturb the system in the same manner. 

2) An exponential noise distribution is assumed, yielding the worst case 
condition of 72 missed and 75 false pulses per hour. 

3) By the use of electronic gating, detection of missed and false pulses 
will be possible. 

4) If a missed (or false) pulse occurs on one of the star signals, the 
system will ignore that signal and use the last stored value of error corres- 
ponding to this signal. 

Resolution of a pulse in time (pulse jitter) and hence the ai^ular error 
associated with detection of the signal depends directly on pulse (or slit) width 
and inversely on signal-to-noise ratio. Using the baseline star tracker parameter 
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values given in the STARS Feasibility Study, the la value of track loop angle 
noise (pulse jitter) is 1.15 arcseconds. 

Figures 4-18 thru 4-21 show the STARS system response when the star tracker 
simulation model is added to the control system simulation of Section 4.2. In 
this case, no disturbance torques are acting on the spacecraft and the torquer 
is still ideal. Figure 4-18 shows the spacecraft attitude time history for half 
an orbit period. A peak to peak error of approximately .20 arcsecond is observed 
in all three axes. This is obviously well within the accuracy requirements of 
the STARS system. Figure 4-19 shows the star tracker's estimate of the spacecraft 
pointing error. If the estimate were perfect Figure 4-18 and 4-19 would be 
identical. However, by a comparison, the measurement is seen to be very noisy, 
yet the spacecraft error is very small. The important point demonstrated here 
is that the error in the star tracker's estimate is filtered by the large inertia 
of the vehicle. Figures 4-20 and 4-21 show the body rates. These rates contain 
a small magnitude, high frequency component in aduition to the required perfect 
pointing rates. 

In conclusion, this simulation shows that the present star tracker design 
will provide the required high accuracy performance. 
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FIGURE 4-18. SPACECRAFT ATTITUDE WITH TYPICAL STAR TRACKER NOISE BUT NO 
DISTURBANCE TORQUES 
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ATTtlUOt tflHOR. ARCStCONDS 
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FIGURE 4-19. STAR TRACKER ESTIMATE OF POINTING ERROR, ZERO DISTURBANCE 
TORQUES 
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FIGURE 4-20. ROLL-YAW BODY 
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PITCH BODY RATE. ARCSECONOS/SECOND 



TIME. THOUSANDS OF SECONDS 


FIGURE 4-21. PITCH BODY RATE WITH TYPICAL STAR TRACKER NOISE BUT NO 
DISTURBANCE TORQUES 


61 


40341-27 (U) 





4.4 STARS System Model 


A signal flow diagram of the complete STARS system, illustrating the 
functional interfaces between the computer, telescope assemb .y, and the 
spacecraft attitude control system is shown in Figure 4-22. The functional 
description of the individual components of the system have been discussed 
previously. 

The total system was simulated on a GE635 digital computer using the 
MIMIC simulation language. The computer simulation flow diagram is given in 
Figure 4-23. The complete STARS system simulation is a combination of the 
simulation models of the precision gimbals, the star tracker, and the control 
system discussed previously, with some additional modelling added for the 
on-board computer requirements. All major disturbances were included. The 
computer listing can be found in Appendix 5.2. 

Figures 4-24 through 4-27 show the complete STARS system response for 
one orbit period. The spacecraft attitude is shown in Figure 4-24 * The 
response consists of an error component at orbit frequency due to the disturb- 
ance torques and a high frequency error component due to the star trackers 
estimate of spacecraft attitude. The maximum peak to peak error in space- 
craft pointing is observed to be well less than half of an arcsecond. The 
star tracker's estimate of the spacecraft attitude is shown in Figure 4-25. 

By comparing this figure to 4-25, the filtering effect of the large mass 
of the spacecraft is observed once again. Figures 4-26 and 4-27 show the 
spacecraft body rates. These rates contain a small magnitude, high frequency 
component, due to STARS sensor noise, in addition to the required perfect 
pointing rates. This simulation shows convincingly that the STARS system 
has the ability of high precision spacecraft pointing. 
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FIGURE 4-22. STARS SIGNAL FLOW DIAGRAM 
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FIGURE 4-23. COMPUTER FLOW DIAGRAM 

















SPACECRAFT ATTITUDE. ARCSECONOS 



0 0.6S 1.30 1.95 2.6 j.* 3.9 4.55 5.2 5.86 6.5 

TIM5, THOUSANDS Lr SECONDS 


FIGURE 4-24. GPACECRAFT ATTITUDE, DISTURBANCE TORQUES AND STAR TRACKER 
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Three addltloaal runs were made using the complete STARS system simulation 
to investigate the effect of gimbal axis misalignment, the effect of a non- ideal 
torquer, ai^ the initial acquisition ability of STARS. 

1) Gimbal axis misalignment 

The effect of an error in the built-in angle between the pitch 
and polar axis of the telescope assembly was considered. A 
three arcsecond error in this angle was assumed to be representa- 
tive of this type of an error. The total STARS system simulation 
described above was run with this three arcsecond gimbal 
misalignment. The system response to this error was not visible, 
thus yielding a conclusion that this type of error will not affect 
the STARS system performance. 

2) Non- ideal Torquer 

A one second deadband region was added to the torquer model shown 
in Figure 4-9. The system response to this error is very small. 

Figure 4-28 shows the spacecraft attitude for this simulation run. 

When comparing this figure with Figure 4-24, an increase in the 
peak to peak magnitude of the pointing error of . 1 arcsecond is 
observea This type of error has a greater effect on spacecraft 
pointing than the gimbal misalignment error, however even with 
this eiror, the STARS control system is well within its accuracy 
requirements . 

3) Initial Acquisition 

The acquisition problem of the STARS system begins when the star's 
image is first seen in the telescope. The image will first appear 
when the error between the telescope boresight axis and the line- 
of-sight to the star is less than .5 degrees (1800 arcseconds). 
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This condition was simulated by giving the spacecraft an appropriate initial 
attitude so that the error in the telescope frame was approximately .5 degrees. 
Tlie STARS system response for this sioiulation is shown in Figures 4-29 through 
4-31. Figure 4-29 shows that the spacecraft pointing error is brought to zero 
in approximately 300 seconds. Figures 4-30 and 4-31 show Che spac .craft's body 
rates, which are also brought to their respective perfect pointing values in 
approximately 300 seconds. 
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FIGURE 4-28. SPACECRAFT ATTITUDE USING CONTROL TORQUER WITH DEADBAND 
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FIGURE 4-29. SPACECRAFT ATTITUDE, INITIAL ACQUISITION 
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FIGURE 4-30, ROLL- YAW BODY RATES, INITIAL ACQUISITION 
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FIGURE 4-31. PITCH BODY RATE, INITIAL ACQUISITION 
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4.5 Conclusions 


A number of conclusions may be reached, based on the analyses, experi- 
ments, and simulations performed in this study. These cover star tracker 
characteristics and performance, the capabilities of the precision gimbal 
system, control loop design, acquisition perfomance, pointing precision and 
pointing stability. 

Star tracker analyses, based on experimental data from several photomultiplier 
tubes. Indicate that without using elaborate electronic processing, the false 
alarm and missed pulse rate to be anticipated lies between 75 per hour and 4 
per hour. Even though these rates appear quite acceptable, an additional order 
of magnitude improvement is available if needed by applying time gating of the 
star pulses as part of the signal processing. 

detailed simulation of STARS gimbal characteristics shows that the gim- 
bals may be positioned with a precision of better than two arcseconds, using 
shaft angle transducers with a resolution of 19 bits. Since the gimbal positioning 
precision appears to be limited chiefly by the resolution of the feedback trans- 
ducer, the prospects of higher pecision in this area (if desired) are very good. 

Simulation of the spacecraft control system shows that a relatively uncom- 
plicated control system design can reduce the effect of various external distur- 
bance torques on pointing precision to negligible values while simultaneously 
keeping pointing jitter due tosensor noise at acceptable levels. It was found 
that a pointing error of 0.5 arcsecond due to disturbance torques and less than 
0.2 arcsecond jitter due to sensor noise is achievable without placing significant 
demands on system gains or bandwidth. At the same time the transient response of 
the overall system appears quite adquate, with acquisition from the edge of the 
field of view occurring in less than 300 seconds: 

The overall conclusion that may be drawn from the results of the STARS 
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System Application Study Is that a relatively straightforward precision pointing 
system design utilizing STARS on an EOS type spacecraft, and considering the 
characteristics of real sensing control hardware, can achieve the desired 0.001 
degree earth pointing capability . The STARS approach should therefore be a 
strong contender In any application where a large payload must be pointed earth- 
ward with high precision from any arbitrary orbit. 
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5. APPENDIX 
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5. 1 Spacecraft Egiiations of Motion 


The spacecraft angular momentum is 


H - I(u 


or, In expanded form 


^RP ^RY 
^PR ^P ^PY 

^YR ^YP ^Y - 


For the EOS configuration specified by GSFC, the only non-zero crossproduct 
of inertia is I (= I™)* Therefore 

Rl iK 


* 'RY^ 




'r'% 'ry% 


,^Y^ ^RY“^ 


The moment equation is 


M = H + (u X H 


OU X H 




U)r U)p 

> »P »Y 

-\'*>R ^ «R“V 
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Thus the expanded nioinent equations are 


^ ^RY “'P % '*' ‘“y ' 

^ ■ ^Y “>Y ^R '*>R "V 

My - Iy Uiy + Ij^Y ^P ’ ^R ‘*’r ‘"p ■ 


ip tUp I^Y 

^RY “V 
^RY "'p "V 
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5.2 COMPUTER PROGRAM LISTINGS 
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% 

IDENT 366»54790»1087D»78294.RnY 

i 

LIMITS 30 

♦ 30X, » » 1 OX ♦ » 1 

$ 

PROGRAM RLH5.NDUMP 

X 

LIMITS 30 

• 30000* flOOOO 

i 

PERM H* 

•MIMIC-SYSTEM 

S 

TAPE 09 

• A25 

CON(Kl.KMl.KIl.KDl) 
CON(TAul»TAUll.Rl.Jl ) 
CON(THETDO*ThETAO) 
DBM(TRUE) 


DTMIN 

1/(2**13) 


dtmax 

DTMIN 


DT 

.01 

F^^ (T.20) 


Sm'1 

F 6*1 (T. FALSE .TRUE* FALSE) 

6m 1 

KROA 

3600*180/3. 14159 

6W1 

KS2 

( K.01*Xl*tC.Ml ) /R1 

6m 1 

OMEG 

( 13.737*360)/24 


SW3 

F SW( REM ( T.,0015 ) .false .TRUE. FALSE) 

Si«y3 

ALPHA 

OMEG*T 

SW3 

THECOM 

ALPHA 

OM3 

thetac 

INGIKI 1*THEC0M) 


IMDTER 

RNG(0.1.75».l) 


BEARER 

RNG(0.1.137..2 ) 


EIIO 

thetac-thetb 


Ell 

K52*EI 10 


P21 

INT(EIl.O) 


Pll 

INT ( ( 10*TAU1+TAU11)*E: 1-F01+P21 .C) 


EOl 

( 10*TAU1»TAU11*EI1+P11)/TAU1 


on 

tE01+Fl)/Jl 


FI 

FSW( ThETAD*. 026.0. -.026 ) 


THETAD 

INT (OII.ThETDO) 


THETA 

INT (THETAD. THETAO) 


THETAA 

krda*theta+oearer 


THETB 

ing(kii*(thetaa+indter; ) 


THEDEL 

thecom-thetaa 


R 

( 1/T)*INT(THEDEL*0) 


RM3 

SQR ( ( 1/T)*INT(THEDEl*T) EDEL.O ) -R*R ) 
PLO( T.RMS) 

SCL(20.0.5.-5) 

END 

200 

.376 

.404 6.25E-3 

.01 

1 

13 1.5 

0 

0 


s 

END JOB 

PRECEDING PAGE BLANK NOT FILMED 

PITCH GIMBAL SIMULATION 

83 



s 

IDENT 

366i 54790 t 10870. 7829A. ROY 

s 

limits 

10.30K..10K..I 

s 

PROGRAM 

RLHS.NDUMP 

i 

LIMITS 

10.30000..10000 

% 

PERM 

H*. MIMIC-SYSTEM 

s 

tape 

09.A2S 

CON(tCF2) 
CON(PHIDO.PHIO) 
PAR(K2»KM2»KI2.KD2) 
PAR(TaU2*TAU12.R2.J2 ) 
PAR (count. SUM) 

DBM (TRUE) 


otmin 

l/(2**5) 


dtmax 

otmin 


DT 

1 

FIN(T.IOOO) 


Sm'I 

FSW(T. FALSE. TRUE. FALSE) 

s#ll 

KS 

{ KD2*K.2«K.M2 )/R2 


S*3 

FSW(REM{ T.90) .FALSE . TRUE .FALSE ) 

S'* 1 

KRDA 

180*3600/3* 14159 

S'*l 

OMEGA 

360/(365.26*24) 

SW3 

BETA 

OMEGA*T 

S*3 

PHICOM BETA 

o* 3 

PHIC 

ING(ia2*PHIC0M) 


BEARER RNG(0»1.137».5 ) 
li^OTER RMG(0*1.75».'') 
******** Elastic stiction model 


OM 3 

CO 

PHI 


S*4 

FSW(ABS( PH ID) -.000001 .FALSE . TRUE . FALSE ) 


SW5 

C0M(SW4) 

S*4 

F2 

n 

S*5 

Cl 

PHI -CO 

S*5 

»*»»*** 

F2 

-KF2*C1 


EI20 

PHIC-PHB 


EI2 

KS*EI20 


P22 

INT(EI2.0) 


P12 

INT( ( 10*TAU2+TAU12)*EI2-E02+P22.0) 


E02 

( 10*TAU2*TAU12*EI2+P12 ) /TAU2 


012 

{ E02+F2)/J2 


?HID 

INT (CI2.PHID0) 


PHI 

IN"^ (PHID.PHIO) 


PHA 

KRDA*PHI+BEARER 


°HB 

ING(ICI2*{PHA+INDTER) ) 


PHIDEL 

PHICOM-PHA 


R1 

( 1/T)*INT(PHIDEL.0) 


RMSl 

SCR( ( 1/T)*INT(PHIDEL*PHIDEL.0) -Rl*Rl ) 


print 

FSW(REM( T.IO ) .FALSE. TRUE. FALSE ) 

PRINT 


OUT (T. PHA. PHICOV. PHI del. RMS l ) 
PLO(T. PHA. PHICOM) 

SCL ( 1000.0.50.0) 

PLOIT. PHIDEL) 
SCL'IOOO.O.IO.-IO) 

PL0(T.RMS1 ) 

SCLdOOO. 0.5.0) 

END 

182.63 


POLAR GIMBAL SIMULATION 

0 

0 
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1.5 .376 .404 

.01 1 13 

1 0 
$ END JOB 


6.25E-3 

.05 
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STAI8 systm wmL 



% 

IOENT 36604750 ♦10H7D*7S294*r0Y 



# 

L1HII3 6a»3QK»*10X**l 



» 

PRCORA ‘‘i i\i.H3 * f4Dy i*.P 




klK|l3 65*300uQ» 




PcR*h n«*;.llNUC-SYSTi,M 



• 

tape 09.A23 



H ali-ULAR TSAC^iNs* aTTITiOC kEFcRCNCL <;Y5TiM 


STARS — 


OATA Input 




constants 



KC 

u I stance from vehicle to Earth 


‘ iLuS 


MAONifUOt Of lARTHS velocity ABOUT 


.<ILES/SEC 

C IO •'*£OLc.CT STi.wLAR A6E!\?.ATI 3ET R0=V0=0) 



OAi'iMA 

0P6IT inclination ANGLE 


OEuFEES 


CC.NCRO»VCtGAs.-*AJ 



ul/.UnL 

V.wii3 1 /ti'i I 3 



K4 

PlTCrt wRlV^R LAIN 


vclt_/vclt 


PilLn '-OlvN lain 


FT-Lu/A:- "5 

l%4 4 

PIiCji IrtiLUwTwLY*^ wmIN 


31 T! /ARCSEC 


Pi 1 1 ^ * w i' S I U.M Ur' I * « 


VOLTl/uIT 

Kr 1 

v»wLU*'*m F R I ^ I 4 C ^ I T C**J AX 4 o ) 


ft-l-s 

i Mw X 

PI TO'. TI-'£ CON^TA.'iT 1 


L L. 4. . 

♦««ii 

=IlLn i I"t CCNbTANl 2 


SECONDS 

Rl 

PiT.-ri XiNCiNL resistance 


CH':£ 

^4 

? 1 TCn • I A 


SL4.G-FT.-ri 


CU.NC Nl *kI 1 •<Ul»<r I J 




lQNCTAuI *TA- ll#RitJl J 




i'wLAR t.-.*. IVlR lai,\ 


/ i- T w / V w L c 

N>‘>& 



F T “L w / r* ‘’*K J 


r'^L^A 1 1 V ww L t Q 3 Y « 4 


w i T 


PwL"i^ L/ A CwN^^'t !%S I Li'« 


y^LTt/blT 

•sr 4 

uLA.STIC FRICTiSii COcFr IeAe-NT (PGLAR 

AX I S } 


i Awi 

P0L*»A T l.-lt CCiNeiAnT 1 


w4» Wwl4ww 


POLAR Tl:.;£ CONSTANT 2 


■jcCcr*w^ 


POLAR 4 1 NO Ins PIESISTANCL 


GHyb 

J2 

POLA.4 axis inlrtia 


T ^ t 


COf4tA2».<’N2*Nl2»XUi*<F2 J 




CON C T AU2 »TA«12 • P2 *J2 C 



C iA3 L 

•">.0.R CC.-.Pw, iSAT I C.N C ' .'.S 1 At»T3 



iirt.L T rtc. 

PITs.!! SI.4 ShL i I as ^RRGR 


^EG?.£LS 

»»t.L.PriI 

POL"R SI”3AL cIAS -RkOR 


4*C.U;;3tlS 

WtoLOA.’. 

GRPlT 1 nClINATIS>« uImS l.^RC.R 


DEURc.ES 

DAJ 

ROLL ( TllESCJPE 1 r- 1 AS FC.R J~STA.R 


uSGREES 

uxa 

ROLL* TEllSCDPL )uIA3 .‘-OP :<-STAP. 


DEC-PELS 

3t0 

YA.<tTcLE3C0Pt)ElAS FOR J-STAR 


DEGP.Ec.3 

oYK 

YAaCTElESLGPcIcIAS rO.R K-STAR 


DEC-PEES 

oZj 

PITCHCTELESCCPElrlAS "CC J-STAR 


OEGi'EES 

UZX 

PITCnCTtLESCCPD&IAS FOR <-5TAP 


DEGREES 


CG.NCOElTHI tlELPnl tOELCA'U 




CON t BX Jf CV J tSZ J *BXE *S YN *£Z<) 



lirti: TO ANoLE C0.'4VERSICN CONSTANTS 



K.i 

SCA.N circle RADIUS 


DEGREES 

P 

SCAN74ING PERIOD 


SEC 

KPiCn 

f'ilSScU PLUw FALSCi PULSc. RATE 


PER HuUR 


CC.N C K.R tP ♦.•.?«« ) 



CuATKOL 3yST£i-i COisSTAiNTS 



A2Y 

CS.NTRQL SYSTt.-i GAIN — YAm 


f T— Lc.3/.\/^g 

K2a 

Control syste'-* gain — roll 


FT-LLS/P.AD 

K<P 

CONTROi. syste.-; GAIN — PITCH 


FT-LE-/RAD 


gg RSPB(X>IK^IUIY OF 1H1 


PAGE m mm 



TAUIY 

lead time Constant — yaw 

SECONDS 

tauir 

LEAD time constant — ROLL 

SECONDS 

TAUIP 

LEAD TIME constant — PITCH 

• SECONDS 

TAu2Y 

LAG time constant — YAW 

SECONDS 

TAuiR 

LAG Tiftt CONSTANT — ROLL 

SECONDS 

TAuZP 

LAG time constant — PITCH 

SECONDS 


C0NlAZY,TAwi*tTAuZY,K2R*TAuXR.TA02R> 
CON IIC2P * TAUlf> * T AU2P ) 

dPACtCRAFT OYNA»nlCS CONSTANTS 


IR 

SPACECRAFT ROLL INERTIA 

slug-ft«ft 

IP 

SPACECRAFT PITCH INERTIA 

SLUG— {"T^FT 

lY 

SPACECRAFT YAn INERTIA 

slug-ft*ft 

ISY 

spacecraft roll-yaw cross inertia 

CONUR>IPtIY»IRY} 

SLUG— FT*FT 

INITIAL 

CONOITIONS 


ALP H AO 

alpha ic 

DEGREES 

&ETAO 

BETA 1C 

DEGREES 

THtT AO 

pitch GIMBAl initial POSITION 

ARCSECS 

THs.Ti/0 

PITCH GIMBAl INITIAL RATE 

ARCSEC/SEC 

PHiO 

polar GIMBAL INITIAL POSITION 

ARCSECS 

PHIOO 

POLAR GIMBAL INITIAL Ri^TE 
CCN(ALPHAO*BETAO) 
CONJTH£TDO»ThETAO*PHlDO»=HlO) 

ARCSEC/SEC 

YAwOlC 

yawdot ic 

ARCSEC/SEC 

ROkOIC 

ROLLDOT IC 

ARCSEC/SEC 

P ITdIC 

PiTCrOOT IC 

ARCSEC/SEC 

YA«0 

YAW IC 

ARCSEC 

RwLO 

ROL IC 

ARCSEC 

PITO 

PIT 1C 

ARCSEC 


PAR ( Y AWO I C f ROLD I C #P I TD I C • YAWO * ROLL 

.PITO) 

TO 

simuuATion initial time 

SECONDS 

IF IN 

SIMULATION FINISH TIME 

ciECONDS 

N 

STEP SIZE s 


PRINT 

printing interval 

SECONDS 


PAR ( T 0 1 TF 1 N »N .PR I NT ) 


STAR selection 


DELTAJ 

DECLINATION OF J-STAR 

OEGREc-S 

deltaic 

i^ECLlNATlON OF K.-STAR 

DEGREES 

LAHDAJ 

RIGHT ASCENSION OF J-STAR 

DEGREES 

lamoak 

RIGHT ASCENSION OF K-STAR 

PAR ( OELTAJ .LAMOAJ »uELT AX .LAMDAX i 

degrees 

lTAK 

RIG.-lT ASCENSION OF X-VELcaCOPE oORcSIGKT 

DEGREES 

ZcTa< 

^EClIhATION of X-TELcSCOPw BCRcSlGHT 

uECRLES 

eta J 

RIGHT ASCENSION OF J-TELEdCO?E BORESIGHT 

DEGREES 


DECLliHATUN OF J-TELESCOPt BORESIGHT DEGREES 

PAR(ZETAJ»ETAJ*ZETA.<fETAKJ 
preliminary calculations 

OBM(TRUE) 



SWI 

FSwtT .FALSE. TRUE. FALSE) 

SWI 

PIE 

3.1AI59 

SNI 

RTD 

180/PIE 

SWI 

C 

186000 

Swx 

krda 

3600*RT0 

SWI 

KA 

kr*krda 

SWI 

KP 

tXDl*Xl*Kf41)/Rl 

SWI 

KS 

CXD2*X2*XM.2I/R2 

s*l 

OMEG 

t U.737*360)/<2A*3600) 

SWI 

OMEGA 

360/1 365.26*2A*3600) 
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alphad 

3^10*3800 


betau 

OMEGA* 3600 


ALPHAB 

OMEO 

i>4yi 

SO«EQA 

SIN40HEGA*T1/RTD) 


CCMEOA 

CL''e?cGA*Tl/!tTO» 

dN| 

COSE 

COSi:3.5/RTOl 


PSi 

^AM 1A 


CPS I 

CGo(PSl/RTO) 


SPG I 

olN(PSI/RTOl 


CGAK?4A 

CG S« GAMMA /R TO 1 

^Ml 

soahka 

SINIGAMMA/RTD) 

;>*1 

co&lto 

COStuELTAj/RTSI 

Sal 

SDELl J 

S1N(DELTAJ/RTD) 

iWtl 

COEL'R 

COSIDELTAiC/RTOJ 

»»i 

SDcLIK 

SlNiOELTAK/RTO) 


CL^<400 

COo(LAMOAJ/kTDI 

4«1 

sla^.oj 

o 1 N 1 L A<SO A J / K T •> } 

i«.l 

CL^’40K 

CCS ( LAi4DAL/>\Tu ) 

^•l 

SLA>iOA 

olNCLAMOAK/RTw) 


CioTAiC 

COSCZETAK/RTL) 


0*CT Ak 

OlNiZtTAX/RTD) 


a.i.1 AJ 

C0o(ZETAJ/RTQ) 

a»al 

OAtTAJ 

SI.N(ZETAJ/kTD) 

Sal 

CtTAK 

COStETAK/RTDJ 

;»ai 

SET AK 

oIN(ETAK/RTDJ 

>*l 

CETAJ 

CCStETAJ/RTDl 

SnX 

setaj 

SIN(ETAJ/kTS) 

s»l 

11 

I Y/ ARC A 

SnI 

12 

IR/ARDA 


13 

IP/ARDA 

swi 

112 

IRY/XRDA 

Kw>\ CONTROL LOGIC 


0*1 

OTMIi^ 

1/(2**N) 

d*l 

UTMAX 

OTMIN 

SWl 

DT 

1 


TI 

T+TO 

FIN(TI *TFIN) 

?CMORY 

fTlHlNGvAHD 

control 


Sn2 

F SA ( RuM (Ti*»2?)fr i-.L SE * TR'Jt. ♦FALSE ) 


S»3 

F Sa C REM ( T I f 90 1 ♦ FALSE * T RUE * FALSE } 


CALPHA 

COo(ALPHA/RTD) 


salfha 

SINCALPHA/RTD) 


ALPHA 

omeg*ti+alphao 

i»*3 

beta 

OHcuA* T 1 +BE T AO 

0*3 

cbeta 

COStBETA/RTDl 

0*3 

SBETA 

SINtPETA/RTD) 

GIMBAL 

ANGLE CC.S.MAND GENERATOR 


thecom 

C ALPHA+DELT he J *3600 


STHETA 

SlNCThETA/RTDl 


CThETA 

cosctheta/rtdj 

0*3 

PH 1 COM 

{ BETA+OELPhI )*3600 

o»'3 

SPHI 

SIN(PHI/RTD) 

S*3 

CPHI 

COotPHl/RTwJ 

hitch 

GIMBAL EGwATlONS 


INOTEP 

RNG{Qtl.73t.l) 


bearep 

RNO(0>l»137t.2} 


thetac 

1NG(XU*THEC0MJ 


Elio 

thetac-thetb 
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£11 

KP*EI10 


P2l 

INT(E11»0) 


Pll 

INTm0*TAUl+TAUllJ*EIl-E01+P21#0) 


EOl 

aO*TAUl«TAUll*EIl+Pll J/TAUl 


Oil 

{E01+F1)/J1 


FI 

F3*{THETAD*XFl*0*-iCFn 


THETAO 

INT(011.THETD0> 


THLTA 

INTITHcTAO.TriETAO) 


T HET AA 

XRDA*THETA^BEAREP 


THtTfa 

INGIiai*CTHETAA+lNOTEPn 

POLAR 

GIM8AL EOUATlQNd 

** elastic friction model 

S'^3 

CO 

PHI 


3*5 

FSW ( ABS (PH ID I-. 000001 .FALSE .TRUE .FALSE > 


SM6 

C0M(S*5) 

S«5 

F2 

n 


Cl 

PHI -CO 

intb 

F2 

-XF2»C1 

»»i;f 

bearer 

RNG(0.1.I37.,S} 


INOTER 

RNG(0.1.7S..-») 

S*3 

PHiC 

ING(Xl2*FHIC0:-n 


EI20 

PHIC-PHB 


EI2 

X3*EI20 


P22 

INTIEI2.0) 


P12 

lNTniO*TAU2+TAUl2J *EI2-C02+P22.0J 


E02 

{ 10*TAU2i^TAw12*EI2+P12 } /TAU2 


012 

(E02+F2I/J2 


PrtiD 

INT(0I2.PH1&0) 


PHi 

INTCPHIO.PHIOI 


PHA 

XRDA»PHI+o£AR£R 


PMB 

ING(<I2*(PHA-HNDT£R) ) 

OPTICAL STAR SENSORS 

S«t2 

ABll 

CThETA»CPHI-STH£TA*CPSI*SPHI 

S«<2 

Ad12 

-STrtETA*CPHl-CTHETA»CPSI*SPHl 

imz 

AB13 

SPHl*SP5l 


AB21 

CTHETA*3PH I+STH£TA*CPS 1 »CPHI 

S^2 

AB22 

-STHETA*SPHi+CTH£TA*CPSI»CPHl 

Sm2 

AB23 

-SPSI*CPHI 

S*2 

AB31 

SThETA*3PSI 

Sn2 

AB32 

CTHETA*5PSI 

S»2 

AB33 

CPSI 

342 

EBX 

YAW+PI R*ROL-ROR*P IT 

3M2 

E3Y 

-P I S»Y AW^ROL+Y AR*P I T 

Sm2 

EBZ 

ROR*YAW-YAR»ROL+P IT 

3t*i2 

RX 

AB11*EBX+AB12»£BY+A813*EB2 

Sw2 

RY 

Ab21*EBX+AB22*E6Y+AB23*EBZ 

Sw2 

RZ 

A83X*EBX+AB32#EBY+AB33*EEZ 

SM2 

RHHXJl 

SZETAJ*C£TAJ*RX+SZETAJ*SETAJ*RY-C* AJ*RZ 

SW2 

RWYJl 

-SETAJ*RX+CETAJ*RY 

Sh2 

RTWXKl 

SZETAX*CETAK*RX+SZ£TAX*££TAk*RY-C2ETAK»RZ 

SH2 

RWYKl 

-SETAX*RX+CETAK*RY 

SM2 

RMiXJ 

LIHiRHHXjl. -1800.18001 

Sw2 

RHHYJ 

LI MIRHHYJI. -1800.1800) 

Sn2 

R^XX 

L1M(RHHXK1. -1800. 1800 


RHHYK 

LlM(RHHYKl.-i8C0»1300) 

PwLSE 

INTERVAL TIMING 

S42 

TXJ 

-P*ASN(RHHXJ/XA)/I2*PIE) 
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T¥J 

P*ASN( RHHY J/KA ) / ( 2*P I E ) 



TK(C 

-P*ASN 1 RHHX X/ XA 1 / 1 2 *P I E i 



TYK 

P*AS^4tRHHY</XA)/(2*PIE^ 

T IMli 

TO 

ANOLfc CONVERSION 

& i 


KPN 

ING(KPRH*tTFIN-T0l/3600> 



MP^ 

I.NGt (TFIN-TCJ*kN0(0*l»l»/P) 



KP2 

ING< (TrlN-TOI*RNOlO»l.2J/P) 



KPS 

ING< (TFI,\-TO)«RNOtO»lt3J/P) 



KPA 

INGt (TFIN-Ta> ♦RNU( 0.1*4 )/PJ 



Ytol 

F £a ( MP 1-KPN » TRUE .TR’JE .FALSE ) 

<>•*2 


YES2 

( HP2^KPN » 1 .RUE .TRUE .FALSE J 

^*>2 


YLS3 

F { K? 3-. .PN * TRUE » T.RuE ♦ FALSE ) 

^mi 


Y£b4 

F Bw { KP4-.'.PN .TRUE * T.R JE .FALSE ) 

.^42 


NOi 

CCKtYESl) 

wa2 


M02 

C0M(Y£S2» 

^«i2 


NC3 

CCVIYES’J 

&*2 


NQ4 

CCM(YES4) 

^<#2 


OITXJ 

RNGiO. 1.15*2) 

-i^»2 


JITYJ 

RNGlO.l.lS.S) 

d«2 


J I > XK 

RN0(0*1.1S*4) 

j«2 


jityx 

RNGtO.l.lS.S) 

•wl 


rhqx j 

<A«SI.bJ-2*PIc»TXJ/P ) +JITXJ 

'Iwl 

.402 

.NC2 


FPl 

rhoyj 

FP2 

K A *5 1 N ( 2 P I w« T Y J / P ) ■*. J I T Y 

f). 0 ^ 

<A»itl;H-2*Piti*TXK/P}+JlTKK 

A.O 

.405 

.'i03 


RHOXik 

FP3 

.'bOA. 


RHC YK 

S A*5 1 N < 2 « p I T V K/ p ) +U I T Y < ^ 

HQA 


FP4 

n.o 

Ve.^>i 


RHOXJ 

TOLmHOXU.DT nIN.TS) 

TfcSi 


FPi 

1.0 

tti»2 


RhOYJ 

TwLlStlOYJ.i^TKlU.TS ) 

I bo2 


FP2 

l.O 

V&a3 


RHOXX 

TUL C RnOXX»OTKIi'i.75 ) 

Y i» j 5 


FP3 

1.0 

¥c.jA 


RriOYX 

tdlcrhgyk.dtmin.ts) 

y£S4 


FP4 

1.0 

GRsIT 

PERruRBATICi 

NS 

d»I 


NUXJ 

OELGAM*BDtLTj*GLAr-;Oj 

tfW 1 


NUYJ 

-D£LGAK*SLAMDJ 

S«4l 


NUXK 

D£lGAK*5DELTK*CLAVDK 

OMll 


.njyk 

-DELGAK*SLA.'.0J 

GTC.LLAR 

aberration 

£i«3 


Vll 

— Alpha b^SALPnA* C CSLT A*COAi‘MA— S oETA^S wAK'-iA ) 



VI 

V ll-ALPnAa*CALP«-'A*( Coc.T«*SwA;';;>tA+SLLT « m>.A ) 

i>a3 


V21 

— ALPHAb*SALPnA*C SoLT A*CGAt«(riA+Cni-T A*3 LAKKA ) 

d«b3 


V2 

V 2 1 "ALPHAS* C Alpha*' ( SBc^T A*SyA^‘..''iA— C l.lt A*C wAi'ii’.A ) 

&«3 


axj 

( VO*COS£* £ O.mEGA+R.S ”■ ( - ' LA.-.D J»V 1 +CL AKO J *V2 ) ) /C 

£W3 


A XX 

' .C*C0SE*S0VEGA+R0»(-SLAyDK*Vl+CLAyDK»V2) )/C 

Sm3 


AYJ 

— ( VO*COHEGA+RG* t SULl.T J*CLA.*'ii> j* V I'+'SDEL T U*LwA.'(i./ j*V2 ) ) / C 

$«3 


AYR 

-tV0*C0MEGA4-R0*(SCELT<*CLAKDX^Vl+SCELTX«SLAK0K*V2 ) )/C 

)£l.£5C0PC BtASES 


d»4l 


XlXJ 

S2ETAJ*C£TAJ*BXJ-*SZETAJ*SETAJ*BYJ-C2ETAJ*BZJ 



XIYJ 

-SETAJ*EXJ+CETA..»SYJ 

Sw| 


XIXX 

S2E)AX«CETAK*BXK+SZETAX*i,FTAK*BYX-C2ETAK*r:'.K 

StMl 


XIYK 

- 3ETAX*dXX+CcTAK*t.YX 

bIaS I 

COMPENSATION 


BA2 


KUXJ 

RHOXJ+ i X IX J+NUX J+AX J ) *3600 
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imZ MUXK RMQXK+JXIXK+NUXK+AXtO*3600 

MUYJ Rri0¥J+«XIYJ+NUYJ+AYj)*3600 

MUYK RHOYK-*-IXIYX+NUYK+AYk)*3600 

tauRDiNATE TRANSFORMATION AND ERROR ESTIMATION 

RU CALPHA*CbETA-SALPHA#Sct"A*COA.MiV,A 

a»2 R12 -SALPHA*CbETA-CAtPHA*SbETA*C4iAMHA 

b«2 R13 SBETA*SGAMMA 

<3h2 R21 CALPHA*sEc:TA-»>SALPHa«C.Sc.TA«CuAMMA 

b«2 R22 -SALPHA*SBETA+CALPHA*CfatTA*CDAK,i.iA 

d*2 R23 -SGAMMA*Ct>cTA 

b«2 R31 SALPHA*SI»AMMA 

i«i2 R32 CAUPHA*SGAM.'.A 

bW2 R33 CGAMMA 

S«i2 RilJ SDcl-TJ*C^.AMDJ^^Kll+SD£LTJ»SEA^iDJ»^<2^-CD£LTJ^^R 

S«i2 R12J SDELTJ*CLA.'’iDj*R12+SD£LTJ*^SLA.-:DJ«R22-CDELTJ*R 

i*2 R13J Sw£LTJ*CLAi';uJ«R13+EOELTJ*SLAMDJ*k23-CuELT J«R 

b*2 R21J -&LAMDJ«k11+CLAMDJ*R21 

b*2 R22J -SLAMDD*R12+CLAMDJ*R22 

Sfc2 R23J -SLAMDJ»R13+CLA.MDJ*R23 

it»2 RllX Su£LTK*CLAi'iDk*Kll+SOELTK*SLA.‘-'.D:«R21-CDLLTR»R 

^*2 R12< SD£LTX*CLANiDK.*R12+5D£LT!(,*SLAMDN->R22-CDELTi<.*-R 

ifeZ R13 k SD£LTK*CLAHDK*Ri3+5D£LTK*SLA>iDN«R23~CD£LTK»R 

R21X -SLAMDA*R11+CLAMD<*R21 

D*2 R22K -SEAMDX*R12+CLAMDX’»R22 

Gw2 R23X -SLAMDK*R13+CLAMDX*R23 

a»»i2 AU KllJ*HllJ+R2lJ*R21J+RIU*RUk+R21k«R2lK 

SK»2 A12 R11J*H12J+R21J*R22J+R11K«R12n+R21K«R22K 

Sj« 2 A13 RllJ*R13J+R2lJ*R23J+RllK*R13N+R21k«R23K 

S42 A21 AI2 

S»2 A22 Ri2J*R12J♦R22J*R22J+^U2K^^R12^+R22^«R22.< 

Si»2 A23 Rl2J*Ri3J+R22J*R23J+R12k*R13K+R22K*R23X 

b*2 A3I A13 

S«2 A32 A23 

S«2 A33 R13J*Rl3J+R23J*R23J+R13.<»R13k+R23K*K23K 

OETAI Aii*(A22»A33-A23»A32 )+A12*CA23*A3i-A21»A3i) 

S»2 OETA 0ETA1+A13*( A21*A32-A22*A31) 

AUIN ( A22#A33-A23*A32)/DETA 

S«(2 A12IN ( A13*A32-A12*A33)/D£TA 

3*2 A13IN ( A12*A23-A13*A22)/DETA 

3*2 A21IN < A23*A31-A21*A33)/DETA 

SW2 A22IN ( A11*A53-A13»A31)/DETA 

0*2 A23IN ( A13*A21-A11*A23)/0ETA 

3*2 A31IN < A21*A32-A22*A31)/DETA 

3*2 A32IN { A12*A31-A11*A32)/DETA 

0*2 A33IN < A11*A22-A12#A21J/DETA 

3*2 ASll A11IN*R11J+A12IN»R12J+A13IN#R13J 

3*2 A312 AliIN*R2lJ+A12IN*R22J+A13IN»R23J 

3*2 AS13 AlllN*RllX+A12IN*R12K+A13IN*Ri3K 

3*2 ASIA AiUN*R2iK+Ai2IN*R22<+A13IN*R23K 

0*2 AS21 A21IN*R11J+A22IN*R12J+A23IN*R13j 

3*2 AS22 A21IN*R21J+A22IN*R22J+A23IN*R23J 

3*2 AS23 A2lIN*RllK-i;,22IN*R12K+A23IN*R13iC 

3*2 AS2A A21IN*R21K+A22IN*R22K+A23IN*R23< 

3*2 AS31 A31IN*RUJ+A32IN«R12J+A33IN*R13J 

3*2 A332 A311N*R21J+A32IN*R22J+A33IN*R23J 

3*2 A333 A31IN*RUK+A32IN*R12X+A33IN*K13K 

3*2 A33A A31lN*R2lK+A32lN«R22K+A33IN<fR23K 

3*2 EAXl AS11»MUXJ+A512*MDYD+A313#?'.UXX+A314».MUYK 
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ii) V>> lo O) Iv! 




LAVi 

AG2X»;iOXJ+AS22«NUYJ+AS23!»MUXF. + AS24«:.:uYf; 

«> 

lA21 

AS 3 l^.'.iUX J+AS32«MUY J+A S 3 3 .•'•■UX K+ A$ 34*:-:uY.. 

i hCL 

SYiTE". 

UEADLAND 


ONI 

FSa (ABC (lAXI) -1 •♦false ♦TP'JF* TRUE) 


0N2 

FSW ( ASS ( CAY 1 ) - 1 • ♦ false ♦ TRUE ♦ TRUE ) 


0N3 

FSu(AbS(EA2l)-l * ♦false ♦TRUE •TRUE) 


OFFl 

CO‘-MOM) 


OFF2 

COM ION 2) 

J 

OFF 3 

COM (ON 3) 

QHl 

EAX 

EAXI 


EAY 

EAYi 

Wl>l> 

tA2 

EA21 

JFFi 

EAX 

0 

QFrZ 

EAY 

0 

QFFJ 

EA2 

0 

S-2 

ERX 

-eax/krda 

i*»2 

ERY 

-EAY/KROA 

i»«2 

ERZ 

>EA2/KRDA 

control 

system 



AC2Y 

INT(2*K2Y#TAUlY*ERX-TY-HNT(K2Y*ERXfU) ♦0) 


AC2R 

INT C 2*K2R*TAUlR*ERY-TR-»-lNT<K2R*ERY#0 ) tO) 


AC2P 

lNT(2»K2P*TAUlP*ERZ-TP+lNTtK2P*ERZ*0)^0) 


TY 

Cx2Y»TAU1Y»TAU1Y*ERX+AC2Y)/TAU2Y 


TR 

( .<2R»TAu1R«TAU1R*ERY+AC2R)/TAU2R 


TP 

( k2P»TAu1Pk-Tau1P*ERZ+AC2P)/TAU2P 

wl jluKcAi'iCc TOR 

UuEo — GRAVITY '■•RADI ENT 1-SOLAR PRESSURE 


T GK 

- 3 . 0» C 0M':G/ R TD ) «■ t OMEG /R TD)»nP-lY)*SIN(?. C*ROP ) 


TGP 

-R.O^CO‘’EG/RTO)*('RMPG/RTO)*{ 1R-IP)*SIN( ?.0*PTR) 


TGY 

o.n 


TcR 

.onooi*r''£{nMPG*Ti/i5TO) 


T SY 

.00CC1«^:;T^'{ OMEG*T ! /RT P ) 


TSP 

^•0 


T A u I S 

TGY-fTSY 


TY'JIS 

TGR-t-TSP 


TZOIS 

tg^+tsp 

^ R 1 

AFT DY.'iA 

•lies — EOb 


DTY 

( { 13-12) I T»'VROL- T 1 2*'yVftWiH;P I T ) /KPDA 


OTR 

( 1 I l-I"* )*»> YXVf*V'OiT+T 1 2*V,’«f'l.*UPTT) AY®PX 


f'l T '■■ 
1 - 

( ( I ?- 1 1 )i".-'YAvi*WPOL+l 1'?*‘(V.YA'-'*WYA'-'-'<'P0L*''RPL ) ) /^®PA 


T1 

TY-DTV+TyP! S 


72 

TR-DTR-fTYDIS 


T3 

TP-CTP+TZDIS 


0-T 

n^i2-n2»iir 


'.*.XD 

( I2»T]-n?»T?) /PET 


.,YD 

( Il<‘T?-n3vTl ) /PET 


.vZD 

T3/I3 

X 

YAV.D0 

.■iXORB + YA'vPIC 

1 

RGLOO 

aypre+rcldic 

5.-.1 

pitdo 

..vZORB+PlTPIC 


>■« YAtS) 

I%T(hXC^YA SO) 


.*^R0L 

INT(WYD^RCLOO) 


•YD it 

IN TCa'ZD ♦PITDO) 


WXORB 

®31i‘0f''FGAi--3600 


WYDRB 

R32^0'’FGA»’600 


W20RB 

R 3 3 »0M E G A 3 6 C 0 +n:-' EG 3 6 0 C 


a'XR 

WX0R5+P I R * .V YCRB-RPR « v;Z CR E 


aYR 

- P I R «■ A X P R B + Y 0 R B + Y A R » A Z 0 R B 


aZR 

RDR » WZ CR B-Y AP » WYCR3-HyZ OR B 


n 



»-* o o 


If. YiR { i.YAW-wXR )+YAR*ROR* I WROL-WYR ) +ROR# { WPI T-WZR ) 

WSBR { WROL-'aYR } - ( WPI T-WZR 1 *YAR 

^fPSR { wROL-»^YR)*YAR+( WPI T-WZR) 

YAW INT (.iYe,R »YA&0) 

ROL INT(WRBR»ROLO) 

f^lT INT(WP5R*PIT0) 

PIR ^^IT/KRDA 

ROR ROL/KRDA 

YAR YAW/KRCA 

Hi I l*WYAv;+Il2*'.'.'ROL 

H2 I?»U'ROL-^n2»WYAW 

H3 I3»WPIT 

H SQR ) 

OwTPs^T DISPLAY 


THEDEL TMECOM-TWTTAA 
•’’HIDEL PH!CO-''-PHA 


OLTP 
OUTP 
OuT P 
OuTP 
OlTP 
OwTP 
OwT=> 
OLTP 
OJTP 
OLTP 
OuTP 
OwTP 
CLT p 
OuTP 
O.JTP 
OUTP 
OLTP 


R ( l/T)*IMT{THEOfL»0) 

RMS SC?(ABS( n /T)4lMT(THEOEL*THFDEL*n)-P*P) ) 

PI < l/T)#IVT(OHir5rL.O) 

R'-'Sl SCR(Aes( n/T)*lMT(t>HlOEL*PHTOEL»a)-Rl«R1 ) ) 

ERSTRX YAW-EAX 

ERSTRV RCL-FAY 

ERSTRZ DIT-EAZ 

OuTP FS.V(REM{TI tPRI NT ) ,F A L SE * TRUE t FALSE ) 

OUT 

OUT 

OUT(TI ) 

OUT(THETA*AL''HA#ThETAA*THPCO»J».THEDEL*RMS) 
OUT(PHI»BETA.PHA.PHICOM,PHIOEL*P«S! ) 
OUTCRHHXJl »RhHYJl»RHHXKl»RHHYd ) 
OUT(RHHXJ»RHhYJ»RHHXK»RHHYK) 
0UT{TXJ»TYJ*TXK,TYK) 

OUTtFDl *FP?»FP3*FPA) 
OUT(NPl»MP2*MP3»MPAtMPfj ) 

OUT (JITXJ»JITYJ»JITXK.JITYK) 
OuT(RHOXJ»RHOYJ.RHOXK»RHOY<) 

OUT (EAXtEAY»EAZ, TY. tr. TP) 
OUTIERSTRX.ERSTRY.EPSTRZ ) 
OUTITY.TR.TP.TI.TP.t-^) 

OUT(WYBR.WRBR.WPBR) 

0 U T ( W Y A, W . WROL » WP I T . Y A W . R OL » P I T ) 

PLO(TI .EAX.EAY.EAZ) 

SCL(TFIN*TC*5.-5) 


PLC (TI .WYAW.WPCL) 

5Cl{TFIN»T0..10.-.1P) 

PLGfTI.WPiT) 

SCL(TFlN»T0.?10.Zn0) 

PLO(TI .YAW.ROLtPIT) 

SCL(TFlN.Tn.5.-5) 

END 


0 

0 

99 



200 

.376 

.404 

5.25F-3 

.026 

.01 

1 

13 

1.^ 


.8 

.376 

.404 

6.25E-3 

1P2.AS 

.01 

1.0 

13 

.05 



0 n 

0 n 0 0 p 

.25 140 
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G»' O £>' GP" 


6.7 


•01 

•01 

loo 


-e.24 

- 8.24 

% 


67 

67 

1000 

0 

0 

0 

3000 
78.25 
78.25 
END JOB 


6.7 .01 67 

6.7 

1000 0 

0 0 0 

0 0 

0 0 0 

50 

- 56.84 305.78 

- 56.84 305.78 


0 

0 
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